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ABSTRACT 

The projected mass of a gravitational lens inside (circular) apertures can be de- 
rived from the measured shear inside an annulus which is caused by the tidal field 
of the deflecting mass distribution. Here we show that also the multipoles of the two- 
dimensional mass distribution can be derived from the shear in annuli. We derive 
several expressions for these mass multipole moments in terms of the shear, which al- 
low large flexibility in the choice of a radial weight function. In contrast to determining 
multipole moments from weak-lensing mass reconstructions, this approach allows to 
quantify the signal-to-noise ratio of the multipole moments directly from the observed 
galaxy ellipticities, and thus to estimate the significance of the multipole detection. 
Radial weight functions can therefore be chosen such as to optimize the significance 
of the detection given an assumed radial mass profile. Application of our formulae to 
numerically simulated clusters demonstrates that the quadrupole moment of realistic 
cluster models can be detected with high signal-to-noise ratio S/N; in ~ 85 per cent 
of the simulated cluster fields S/N ^ 3. We also show that the shear inside a circular 
annulus determines multipole moments inside and outside the annulus. This is rel- 
evant for clusters whose central region is too bright to allow the observation of the 
shear of background galaxies, or which extend beyond the CCD. We also generalize 
the aperture mass equation to the case of 'radial' weight functions which are constant 
on arbitrarily-shaped curves which are not necessarily self-similar. 



1 INTRODUCTION 

As was shown by Kaiser & Squires (1993), the observable image distortion of high-redshift galaxies by an intervening mass 
concentration acting as a gravitational lens can be used to reconstruct the two-dimensional mass map of the deflector. The 
local mean image distortion provides an estimate for the local tidal gravitational field, from which the underlying surface mass 
density can be obtained. Since then, the original inversion method has been generalized in several ways (see, e.g., Schneider 
1996a for a recent review and references). The mass maps obtained in this way are noisy, and the significance of individual 
features is not easily assessed. Kaiser (1995) proposed to use a measure for the mass convolved with a compensated top-hat 
filter. This so-called ("-statistics, or aperture mass, has the advantage that the statistical properties of this mass estimate can 
be easily obtained from the observations. This statistics was then modified to allow for more general radial weight functions 
(Kaiser et al. 1994), and proposed by Schneider (1996b) as a means to search for 'dark' mass halos on wide-fleld images. 

Recently, Wilson, Cole, & Frenk (1996) have investigated the possibility of studying the asymmetry of galaxy clusters with 
weak-lensing methods. As is well known (e.g., Richstone, Loeb, & Turner 1992; Bartelmann, Ehlers, & Schneider 1993), the 
amount of asymmetry in clusters depends quite strongly on the cosmological model, since it signifles the degree of relaxation, 
and thus age, of the clusters. High-density models predict a much later cluster formation epoch, and thus more asymmetric 
clusters, than low-density models of the universe. Wilson et al. (1996) applied the Kaiser & Squires (1993) cluster mass 
reconstruction algorithm to a sample of numerically generated cluster models and investigated the quadrupole moment of the 
area enclosed by a selected isodensity contour. Doing so, they showed that this measure indeed discriminates well between 
high- and low-f2 universes if applied to a fairly small number of massive clusters. 

Here we derive mass multipole moments in circular apertures in terms of the shear in an annulus. As in the case of aperture 
masses, the chief advantage of aperture measures for the multipole moments is its direct relation to observable quantities. 
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which allows to unambiguously quantify the accuracy, and thus the significance, of the result. Multipole moments can also be 
determined from mass maps derived from cluster mass reconstructions (e.g., Kaiser & Squires 1993; Seitz & Schneider 1995; 
Bartelmann ct al. 1996; Wilson et al. 1996). Apart from being technically more involved, the inference in this case is indirect 
because multipoles are determined from the reconstructed mass maps rather than directly from the shear data, it does not 
allow to give simple error estimates, and it depends on artificial parameters like a smoothing length. 

In contrast to the aperture mass, the weight functions (filters) do not need to be compensated. With two different 
methods, the aperture multipoles in terms of the shear field are derived in Sect. 2, leading to two apparently different 
expressions. It is shown in Appendix A that these two expressions are equivalent, and they can be suitably combined to 
allow for a large flexibility in the choice of the weight function. Several examples for this are presented. Alike the aperture 
mass, the mass multipoles inside a circle can be derived from the shear in an annulus outside that circle. In particular, as 
will be shown explicitly in Appendix B, the shear in a circular annulus is determined by the mass in this annulus, and the 
multipole moments of the mass inside and outside the annulus. We then investigate in Sect. 3 the signal-to-noise statistics of 
the aperture multipoles, using as a simple example a quasi-elliptical mass distribution. In Sect. 4 we present a preliminary 
analysis of the quadrupole moments of numerically generated cluster models to demonstrate the applicability of our method to 
realistic clusters. As discussed in Sect. 5, the application to a larger sample of cluster models drawn from different cosmological 
simulations provides an alternative approach to the discrimination between different cosmological scenarios. 

The methods used for deriving the aperture multipoles are then used in Appendix C to generalize the aperture mass to 
weight functions which are constant on a set of nested arbitrarily-shaped curves which are not necessarily self-similar. 



2 APERTURE MULTIPOLES 
2.1 Definitions 

Consider a mass distribution, e.g., a cluster of galaxies, at redshift «d, with surface mass density E(5). The dimensionless 
surface mass density k{x) is defined as usual in gravitational lensing, 

«(^) = ^ , (1) 
with the critical density 

E 12-) 

and the D's denote the angular-diameter distances between observer, deflector and source, and between deflector and source 
(for notation, see Schneider, Ehlers, & Falco 1992). For simplicity, we will provisionally assume that all faint galaxies are at 
the same redshift z^. The generalization to sources distributed in redshift will be discussed in Sect. 3. 

We define the tensor of quadrupole moments of the two-dimensional mass distribution k{x) with respect to the point xo 

as 

Qij ^ J <fxK{x + xo)w{\x\)xiXj, (3) 

for i,j £ {1, 2}, where w{x) is a radial weight function. Combining the trace-free part of this tensor into a complex number 
Q = Qii — Q22 + 2iQi2, this can be written as 

Q= dxx'^w{x) / dip e^^'^ k{x + Xo) , (4) 
Jo Jo 

whereas the trace of Qij becomes 

f°° 3 ^ ^ 

M:=Qii+Q22= / dxx w{x) / d(pK{x + xo). (5) 
Jo Jo 

These equations are easily generalized for higher multipole orders. We define the complex n*'^-order multipole moment by 

POO /■27r 

Q("^= / dxx"+''w{x) d<pe™^K(f -l-fo) , (6) 
Jo Jo 

and the 'mass moments' 

M*"' := / da;a;"+'w(a;) / d^f ti{x + xo) ■ (7) 
Jo Jo 

Obviously, Q^^' = Q and M^^^ = M. 
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2.2 Multipole moments in terms of the shear 

Let i^{x) denote the deflection potential, which is related to the surface mass density k through the Poisson-like equation 

= = i (V'.ii + V',22) , (8) 

where indices preceded by a comma denote partial derivatives. The two components of the local shear "f{x) = 71 + 172 are 
then given by 

71 = ^ (^,11 — ^,22) , 72 = V',12 • (9) 

It was shown by Kaiser (1995) by suitably combining third-order derivatives of ?/) that the gradient of k can be written in 
terms of derivatives of the shear components as follows: 

. (10) 

\72,i — 71,2 / 

This relation can now be used to express the multipole moments defined above in terms of the shear 7. Since the shear is 
directly observable from the distortion of the images of faint background galaxies (at least in the case of weak lensing, i.e., 
K <C 1), the multipole moments can be written directly in terms of observables. 

For later convenience, we define the tangential shear 7t(a;; xo) o.'nd the radial shear ^r{x; xg) at position x = xe"^ relative 
to position xq by 

7t(f;fo) = - [71 cos(2v5) + 72 sin(2((9)] = -5R [7(5 + fo) e"^"^] , 
7i.(f;fo) = - [72 cos(2v3) - 71 sin(2v3)] = -Si [7(5 + xo) e"^"^] , 

(11) 

where K(z) and Q{z) denote the real and imaginary part of the complex number z. 

2.2.1 First derivation 

We integrate the (^-integral in (^) by parts to obtain 

Q<"' = - / dxx'^+^wix) / dvJe'"*'-^. (12) 
Jo Jo 

The derivative of k with respect to ip can be substituted after transforming ( |lo| ) to polar coordinates. Subsequently, terms 
containing partial derivatives with respect to x and <^ can be partially integrated with respect to x and </?, respectively. This 
leads to 

Q'"' = Q*,"' = / Axx"'^^w{x) I d<^e""^ 7t(a;; ^o) H / da; \jix'^^^w{x) + x'^^'^w' {x)\^ I d(^ e""^ 7r(x; a;o) , (13) 

Jo Jo Jo Jo 

where w' (x) is the derivative of w{x), and w{x) has to be continuous and piecewise differentiable. In order to get rid of the 
boundary terms in the partial integrations leading to (^^, w{x) has to satisfy the conditions 

I7I x"''^^ w{x) for a; ^ and x ^ 00 . (14) 

2.2.2 Second derivation 

Alternatively, we can integrate (o) by parts with respect to x and find 



2-n" c. 

Q<-"^^-l dxxW{x) I dipe""^^, (15) 



/o ^0 
where W{x) is related to the weight function iu(x) by 



xW{x)^ / dyy"+^w(y). (16) 
Jo 

Substituting into ( [l5| ) the radial derivative of k, from (^o|) transformed to polar coordinates, we obtain 

poo p2iT pQC p2n 

(3'"> = Q^"'=/ dx [2W(x)-x''+^w{x)] / d(^e™^7t(f;fo)-i?i / dxW{x) dip e"'"^ -yAx; xo) . (17) 



/o Jo Jo Jo 

The boundary terms in the partial integration leading to ([l5|) vanish if W{x) satisfies 
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xW{x)n-^Q for a:: and x ^ oo . (18) 

We have thus derived two expressions for Q'"' which are of quite different form. We prove their equivalence in Appendix 
A, using the fact that the two components of the shear are not mutually independent, but related via the underlying deflection 
potential. Another way to see the equivalence is to note that the curl of the vector field on the right-hand-side of (|lo|) must 
vanish, which is a condition not used in the derivation of either ( p^ or (p^). 

After a sequence of partial integrations similar to those used before, the mass moments (^) become 



M<") = / Ax[2W{x)-x''+''w(x)\ / dv57t(£;fo), (19) 
Jo J a 

which is the mass aperture equation of Kaiser et al. (1994; see also Schneider 1996b). 



2.3 Aperture multipole measures 

The fact that we have two different, but equivalent, expressions for the multipole moments Q*-"-* provides some freedom in 
choosing the weight functions 10(2:) and W{x) suitably. Recall the situation for the aperture mass measures as discussed in 
Schneider (1996b). There, from shear data in a ring with x £ [uR, -R] with v < 1, the aperture mass measure m = can be 

determined, provided the weight function w^x) satisfies J dx xw(x) = 0. This condition arises because the shear determines 
the surface mass density only up to an additive constant, and the normalization condition on w guarantees that this unknown 
additive constant drops out of the aperture mass. As an aside, we generalize the aperture mass formula to arbitrarily shaped 
apertures in Appendix C. 

Similarly, we assume here that shear data are given in the ring x £ [vR, R] , where in practice the outer radius R is 
determined by the size of the CCD, and an inner radial distance from the cluster center may be required if the cluster 
contains bright central galaxies which outshine the images of faint background galaxies. Then, Q^"' is a local estimator for 
Q("\ in the sense that if the weight function w{x) is nonzero only in a specified radial interval, the estimate Q^"' requires 
shear data only from that interval. However, it should be noted that the boundary terms in the integrations by parts with 
respect to x leading to ( |l3| ) vanish only if w^x) = at the boundaries of the interval. 

In contrast to Q^"\ Qi"^ is a non-local estimate of Q'"\ for it requires shear data from outside the interval, unless W{x) 
also vanishes outside the specified interval. It follows from (|l|) that W(x)^0 outside an interval imposes integral constraints 
on w{x), whereas w{x) is an arbitrary function in the estimator Q^"' (except that it has to vanish at the interval boundaries). 

To simplify notation in the remainder of the paper, we introduce the abbreviations 

gt{x) ~ I dvpe""^ 7t(f; fo) , gi{x)~i I dt/je"'^ 7r(f; £0) • (20) 
Jo Jo 



2.3.1 Local estimate of Q*^"^ 

Let wi (x) be a weight function which vanishes outside the interval x £ [I'R, R\ , and let wi (vR) — wi (R) — 0. Then, the local 
estimate of the multipole moment Q'"' corresponding to this weight function is given by 

qJ"' = / dx x"'^^ wi{x) gt{x) / dx x"'^^ [n'Wi{x) + xw'i{x)'^ 9rix) . (21) 

Ji^R ^ JvR 

This equation is applicable to the case when multipole information is to be determined on the circle x < R (i.e. 1/ = 0) or on 
the ring uR < x < R from shear data measured in the same regions. 



2.3.2 Non-local estimates of Q*^"^ 

We now want to address the following question. Suppose data are available within the radial interval vR < x < R only, can 
we estimate the multipole moments for x < vR and x > R? Since and Qi"' are equivalent estimates of Q'"', one can 
combine the two with weights a and (1 — a), respectively, to obtain 



Q^"^ ^aQ';'> + {l-a)Q, 



(n) 



dx [(2a - l)x"+^ w{x) + 2(1 - a) W{x)] gt{x) 



+ 



da; 



" wix) H X w {x) — n(l — a)W{x) 

n 



(22) 
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In the light of this equation, we ask whether we can choose the weight a and the weight function w{x) such that the integrands 
in (|2^) vanish outside [i^R, R], but with w{x) not identically zero outside that interval. If that were possible, we could deduce 
some information on the multipoles of the deflector in those regions where no shear information is available. 

The requirement that the brackets in the two integrals of ( [22[ ) vanish leads to two solutions for the weight function w{x) 
and the corresponding weights a, 

w+(x) — const, with a+ = —. ^ ; w-(x) oc x~^" with q_ = —. ^ , (23) 

' 2(n + 1) 2(n - 1) ^ ^ 

where the latter is valid only for n > 3. 

We first consider the multipole moment within the circle x < vR, i.e. interior to the region of the data. Thus, we search 
for a weight function 'W2{x) which vanishes for x > R, but which is finite for x < uR. Hence, 'W2{x) must be of the form ( |23| ) 
for X < vR. Since ■W-{x) does not satisfy the conditions ( |l4| ) and ([l^), W2{x) = const, for x < vR, and we choose W2{x) = 1 
for x < vR without loss of generality. In order for W2{x) to vanish for x > R, we further require 

R 

dx x"+^ ■W2ix) ^ . (24) 





Hence, the 'internal' multipole moment corresponding to the weight function W2 with the aforementioned properties 

becomes 



Qr = 2(1 + n) y "^'^ I + 2n)H/2(a;) - 2x"+'w2ix)] g,{x) + [nx"+'w2{x) + x^+^w'^ix) - (2 + n)nW2{x)\ g,{x)^ , (25) 

where we have inserted the weight a+ from (p3|). 

Next, we consider the multipole moment outside the circle x = vR, i.e. exterior to the region of the data. We therefore 
require a weight function which vanishes for x < vR, but is non-zero for x > R. This weight function, called W3(x), needs to 
satisfy the condition - see (|l^) and (^^ - 

dx x"+^ W3{x) ^ , (26) 

'vR 

and if it behaves like w{x) oc x^'^" for x > R, then the multipole moments become, for n > 3, 

Qi"' = I [2a;"+^u-3(x) + (2n - A)W3{x)] g,(x) + [nx^+'wi{x) + x'^+^w'^ - n{n - 2)W3{x)\ . (27) 

This equation is not valid for n — 2 because then w oc x^* for x > R, and W{x) becomes logarithmic. We have not been 
able to find a relation for the external quadrupole moment. We show in Appendix B that the shear in the annulus caused by 
matter not in the annulus is completely determined by the set of multipole moments of the matter inside the annulus (i.e., at 
X < uR) and the multipoles of order 71 > 2 of the matter outside the annulus. 

Of course, the different local and non-local expressions for the mass multipoles can be combined appropriately. As a 
simple example, consider the case that we want to infer Q'"' for < 2; < _R. We then need to combine the internal multipole 
moment (|2^) with the local multipole moment (^l|). Thus suppose we have a weight function w^x) which is constant for 
< a; < uR and smoothly falls off to zero for x —> R. Then, w{x) can be decomposed into 'w{x) — wi{x) + W2{x), where 
wi{x) vanishes for x < vR and x > R, and W2{x) is constant for x < vR, vanishes for x > R, and satisfies the constraint (p^. 
Since the aperture multipoles Q'"^ are linear in the weight function w{x), different expressions for Q'"' can be appropriately 
superposed. To use a weight function of the form just discussed, one would have to insert wi{x) into ( pl[ ) to obtain Q^"-' and 
^2(2;) = w(x) — wi{x) into ( pst ) to obtain Q^^\ and the multipole moment would be the sum of the two, Q("> = Q^"' + 



12 



3 MULTIPOLE SIGNAL-TO-NOISE RATIO 

The aperture multipole moments calculated in the previous section are all of the form 

nR 

Q("'=/ dx[hix)gt{x) + K{x)grix)] . (28) 

JiyR 



With (|20|), we can rewrite (g^) as 

Q("'(£o)= /d^:re"'^ ( ^o) + i^^r{x; xo)] , (29) 



where we have re-introduced the explicit dependence on the center xo of the aperture. 

In order to apply this equation to data, one needs to approximate the integral by a sum over individual image ellipticities. 
Let ei be the complex ellipticity of the i-th galaxy image at position Xi. As in Schneider (1996b), the modulus of the ellipticity 



© 0000 RAS, MNRAS 000, 000-000 



6 Peter Schneider & Matthias Bartelmann 



for an elliptical image is defined as |e| = (f — r)/(l + r), where r is the axis ratio, and the phase is twice the angle enclosed by 
the major axis and the positive xi-axis. In the case of weak lensing, k ^ 1, is an unbiased estimator for the shear '^(xi). 
Since then 7 depends linearly on the distance ratio Dds/Ds, the following relations are valid in the case of weak lensing for 
a redshift distribution of the sources, with Dds/Ds replaced by (Dds/Ds), the average over all galaxies used. Since we will 
mainly be dealing with those regions of the lenses where k ^ 1 (e.g., by choosing the aperture radii appropriately), we write 
the discretized version of as 



JV 

'(Xo) = - > e ' 

j=i 



where n is the number density of galaxy images in the annulus. We have introduced polar coordinates {yi,&i) such that 



(30) 



Xi = Hi 



cos i9i 
sini^i 



+ Xo 



and the tangential and radial components of the image ellipticity are defined in analogy with (uAh as 



eti — — 5R ( e 



eri = — ^ [ei e 



In the absence of a lens, the expectation value of Q'"' is zero, and its dispersion is 



n2 2 ^ 



+ 



where we have used that, in the absence of lensing. 



{eti^ij) = , 

and CTe is the dispersion of the intrinsic source ellipticity. The ensemble average of the dispersion is 



n 



2xk dxk 

(1 - !/2)i?2 



Ay 



We can therefore obtain an approximation for the signal-to-noise ratio of the multipole moments, 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



To obtain an estimate for the expected signal-to-noise ratio of the quadrupole moment of an elliptical mass distribution, 
we consider a singular isothermal ellipsoid, with surface-density distribution 



K,{x) = 



(37) 



2^{l-r,Yxl + {l + r^yxl 

It has an axis ratio of r = (1 — 77)7(1 + 77), and x^ is the characteristic angular scale, which we choose as the Einstein radius 
of the corresponding singular isothermal sphere with velocity dispersion a„, 

\ Ds 

becomes 



For this mass distribution, the quadrupole moment 



(38) 



^(2) ^ XECji) 

2 



da; x^" w{x) 



(39) 



where, for moderately small 77, C{ri) = nrj + 0{rj^). If we choose ( ^ ) for the determination of the quadrupole moment, the 
signal-to-noise ratio (js^) reads 

J dx x^ wlx) 



j(2) ^ C(t?) XeVu 



^ J dxx^^ {2w^ + xww' + x^w'^/4:) 



(40) 



One can now try to find a weight function which maximizes Sc for the mass distribution under consideration. From the 
Cauchy-Schwarz inequality, it is easy to find that the optimized weight function is w{x) oc x~^. If this is inserted into (^o[), 
one sees that the integration has to be terminated at both ends, say at Xmin and Xmax. We then find 



Vln(a;max/a;,„in) ~ 4.5 ( ) ( TtI 



0.2 



30 arcmin 



1/2 



lOOOkms-iy \ 



(41) 
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where is the second step we assumed that (a;max/a;min) ~ 10, and we used the approximation for C{ri) as given above. 

The optimal weight function cannot be simply chosen as a;"'', since we have to satisfy (p^. For the simulations presented 
below we therefore choose a weight function which behaves approximately like (and, for higher multipole moments, like 
^-(n+i)^^ but with appropriate behaviour at 2: = uR = and x = R, 

W,) = I I {n + l)R-(x-R) 

Using this weight function, we have made simulations by assuming a mass distribution of the form 

-3/2 



«(^1 = ^(^!^^„)2 [2 (1 + 0^? + 1-1'] 



2 2 
2 X2 

' (1-77'P (1+77': 



(43) 



which corresponds to an elliptical deflection potential of the form 

^ = (^^' + (T^ + (ifW) ' ' ^''^ 

with axis ratio (1 — f?')/(l + v')- In tl^s limit of vanishing ellipticity, rj' — 0, and zero core size, Xc = 0, (^3|) becomes the 
mass distribution of an isothermal sphere, with x-e given by (^8|). For jxj ^ Xc, the isodensity contours have an axis ratio of 
(1 - r7')V(l + ^7')^ which is 0.74 (0.55, 0.40) for 77' = 0.05 (0.1,0.15). Galaxies with (Dds/A) = 0.7 and a surface density 
of 30arcmin~^ were distributed and sheared according to the mass distribution (^^. The galaxies were assumed to have a 
Gaussian ellipticity distribution with = 0.2. In Fig. |^ we show the result of one such simulation, with = 1000 km s~^, 
a^min = 0.1-R, R — 5', Xc — 1.5a::E (i.e., the lens is subcritical), and four different values of rj' . We have plotted there the 
cumulative probability for the signal-to-noise ratio 

S=^- (45) 

as given by ( |30[ ) and (^). As can be seen from Fig. ^, the typical signal-to-noise ratio of the quadrupole moment for a lens 
with 77' ^ 0.1 is considerably larger than the median of the S-distribution for a spherical lens. That means that for these 
mass distributions the quadrupole moment can be detected at a statistically significant level, whereas for the distribution 
with rj' = 0.05, there is a considerable overlap between the probability distribution of the elliptical lens and the corresponding 
spherical mass distribution. 



4 APPLICATION TO NUMERICAL CLUSTER SIMULATIONS 

Having seen in the previous section that the aperture quadrupole from gravitational shear can significantly be detected for 
moderately elliptical idealized lens mass distributions, we proceed here to study the aperture quadrupole for numerically 
simulated clusters. From a large sample of gas-dynamical cluster simulations performed within the standard CDM cosmogony 
(for a detailed description of the sample see Bartelmann & Steinmetz 1996), we select a subsample of 13 clusters, all at 
redshift Zd ~ 0.3. Since we can project them along three independent spatial directions, we have in total 39 surface-mass 
distributions from these clusters. Assuming sources at redshift Zs — 1, their maximum convergence values range between 
0.27 < Kmax < 1.40, with the median at Kmax — 0.73, and one third of them are critical lenses, i.e., they produce critical 
curves. 

For each of the 39 simulated cluster fields, we simulate 200 source galaxy distributions with a galaxy density of 40 arcmin"'^. 
Since the fields have a side length of 5', the average number of source galaxies per cluster field is 10"^. Their positions are 
distributed randomly, and they are assigned an ellipticity drawn from a Gaussian distribution with dispersion CTe = 0.15. 

We then determine from each simulation the aperture quadrupole estimate Q'^-* according to ( |30| ) and the aperture 
signal-to-noise ratio according to (^), using a weight function w{x) of the form given in (|42[). There, R — 2'.5, and we choose 
Xmin ~ 1'. Since w{x) has the dimension (length)"^, the quadrupole moment has the dimension (length), and we give it in 
units of arc minutes in the following figures. 

Figure |^ shows for each cluster in the sample the fraction of the 200 lensing simulations for which the aperture quadrupole 
is determined with a signal-to-noise ratio (S/N) larger than some threshold {S/N)o, as a function of the modulus of the true 
cluster quadrupole as defined in (^). The difference between the four panels in the figure is the signal-to-noise threshold 
{S/N)o. In panels (a,b,c,d) we have chosen, respectively, {S/N)o ~ {2,3,5, 10}. All panels show that the significance of the 
aperture quadrupole measurement increases with increasing intrinsic quadrupole jQtrueli expected. While the aperture 
quadrupole can be determined with (S/N) > 2 for all clusters with IQtrud <^ 0'.05 (panel a), IQtruc ^ 0^1 is required for 
(S/N) > 5 (panel c), and only a small fraction of clusters allows to determine the aperture quadrupole with (S/N) > 10 
(panel d). For further information, the diff'erent symbols denote the maximum convergence value Kmax of the clusters. We use 
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Figure 1. The cumulative probability distribution P(< S) for a signal-to-noise ratio, as defined in ([45[), being greater than S. For these 
simulations, a galaxy density of 30/arcmin^ has been assumed, with an intrinsic ellipticity dispersion of a,_ = 0.2 and a mean distance 
ratio {D^^/Ds) = 0.7. For the lens, a mass distribution of the form was assumed, with t7„ = 1000 kms^^, core size = 1.5xb, and 
four different ellipticity parameters rj' = (solid curve), 0.05 (dotted curve), 0.1 (short-dashed curve) and 0.15 (long-dashed curve) 

triangles (squares, hexagons) for Kmax < 0.7 (0.7 < Kmax < 0.83, 0.83 < Kmax), where the intervals are chosen such that they 
encompass one third of the cluster sample each. The distribution of these symbols along the curves shows that the aperture 
quadrupole can only be determined significantly if ft;max is not too small. 

While the shear 7 enters into the formula for the aperture multipoles (^^, the distortion of background objects measures 
the reduced shear, which is a combination of shear and convergence, <; = 7(1 — k)~ , rather than the shear. In the case of 
weak lensing, g ~ 7, and then the measured galaxy ellipticities provide an unbiased estimate for 7. If lensing is not weak, 
however, the observed aperture multipole ( |30[ ) is a biased estimate of the true multipole. Since k > 0, Q'"' as determined 
from ( po| ) overestimates the true aperture multipole. This is demonstrated in Fig. ^, where we show the distribution of 
aperture quadrupole measurements relative to the true aperture quadrupole for three clusters with approximately equal true 
quadrupole, but with different maximum convergence values Kmax. 

As a test, we have verified that the systematic bias in the quadrupole measurement disappears when we distort the galaxy 
images with the shear alone rather than with the reduced shear. Apart from the bias, the distribution of points in Fig. ^ 
shows that the scatter in the quadrupole measurements is fairly low. 



5 DISCUSSION 

In this paper we have investigated the possibility to derive multipole moments of the projected mass distribution of clusters 
with weak-lensing techniques. Following the same ideas underlying the aperture mass measures (Kaiser 1995, Kaiser et al. 1994, 
Schneider 1996b), we have derived several expressions for the mass multipoles of clusters in terms of the shear distribution in 
a circle or an annulus around the cluster center. The shear distribution is observationally accessible through the ellipticities 
of faint background galaxies. In contrast to the aperture masses, the weight functions used for the multipoles do not need to 
be of zero total weight, because the additive constant in the mass reconstruction arising from the mass-sheet degeneracy is 
irrelevant for multipoles. 

Our different expressions for the multipole moments can be adapted to a variety of observational situations: 

- Expression (21) permits to determine the multipole in a circle or a ring from data given in the same region. 

- Expression (25) allows to determine the multipole inside a circle from the shear data in an annulus which excludes the 
central region of the circle. This is relevant in cases when the central part of a cluster is dominated by bright cluster galaxies 
which outshine the faint background galaxies. 
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Figure 2. For each cluster in tlie sample with intrinsic quadrupole as defined in (Q), the fraction of the 200 lensing simulations 

is shown for which the aperture quadrupole is determined with signal-to-noise ratio (S/N) larger than some limit {S/N)o, where 
(S/N)o = {2, 3, 5, 10} in panels (a,b,c,d), respectively. The triangles (squares, hexagons) distinguish between clusters with low (medium, 
high) maximum convergence Kmax (see text for details). The curves show that the aperture quadrupole can be determined with (S/N) > 2 
for all clusters with IQ^^^gl ^ 0^05 and with (S/N) > 5 for all clusters with |Qt^^^,| <; 0!l, while only few clusters allow to determine the 
aperture quadrupole with (S/N) > 10. 



- Expression (|2^) yields the multipole in a ring as before, but also outside the region where data are given. Hence, the 
multipole of a cluster can be determined even if the shear data do not cover the entire cluster. 

- As discussed at the end of Sect. 2, these different expressions can be appropriately combined to match specific observational 
circumstances. 

Analytic estimates and numerical experiments using simulated cluster mass distributions demonstrate that the quadrupole 
of massive clusters can be detected in a statistically significant way. As an aside, we have generalized the mass aperture relation 
to apertures of arbitrary shape, both for self-similar curves of constant weight (see also Squires & Kaiser 1996) and for arbitrary 
curves of constant weight. 

The weight function 'w{x) used in our numerical studies might not be the most appropriate choice. In particular, this 
weight function, when applied to the multipole formula (^3|), uses shear information from the entire circle |2;| < i? and thus 
does not punch out a possible region of strong lensing in the clusters. Therefore, the measured image ellipticities are not an 
unbiased estimator of the shear in the inner part of the cluster, which systematically biases the aperture quadrupole moments 
high, as seen in Fig. ^ However, as explained at the end of Sect. 2, one can construct weight functions which avoid this 
problem and which still can be optimized for the detection of multipole moments of quasi-isothermal mass distributions. 

Whereas our numerical results on clusters drawn from gas-dynamical large-scale structure simulations are encouraging 
and give us faith in the usefulness of the aperture multipole measures, they should be considered as a preliminary application 
only. The cosmological application of our new approach to quantifying the asymmetry of cluster mass distributions will require 
the generation of a sample of clusters simulated in different cosmological models. The comparison of the aperture multipole 
moments of clusters formed in different cosmologies will then show whether they are sufficiently powerful indicators of the 
cosmic density parameter fl, and how our new method compares with the one suggested by Wilson et al. (1996). To compare 
clusters drawn from different cosmological simulations, not the quadrupole moments themselves should be compared because 
they depend on both the clusters mass (or velocity dispersion) and the degree of asymmetry, but normalized quadrupole 
moments. As such, the ratio of the quadrupole moment and the aperture mass seems to be an appropriate dimensionless 
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Figure 3. Distribution of aperture quadrupole estimates as obtained from (|30|) witli simulated galaxy cUipticitics (crosses), together with 
the true aperture quadrupole (filled squares), for three cluster models with approximately equal true quadrupole but different maximum 
convergence values Kmax- In panels (a,b,c), Kmax = {0.6, 0.7, 1.0}. The figure shows that the measured aperture quadrupole overestimates 
the true one because image ellipticities measure 7(1 — K,)~^ > 7 rather than 7 only. 



measure of the degree of asphericity of a cluster. A detailed study of these issues, together with a comparison with the Wilson 
et al. method, will be the subject of a forthcoming study. 



APPENDIX A: EQUIVALENCE OF qI^\X) AND Q*f' (X) 

We will now explicitly show the equivalence of ^s]) and (^^, starting from a general multipole expansion of the deflection 
potential. Choosing coordinates such that xq — 0, we can write 

00 

l/jix) = ^ fm{x) COs{mip + ifm) (Al) 
m=0 

the function fm(x) are assumed to have continuous second derivatives. Using (^), (^), and (|ll|), we can then write 71, 72, and 
K, and hence also 7t, and 7r in terms of the functions fm{x), 



= 2 



m=0 



1 / , \( fill \ f'mi^) , 2fm{x) 

- cos{mip + V5m) I fm{x) h m 



= ^ m sm{m(p + (pm) 



.fL{x) fm{x) 



(A2) 



Using further the relations 

2tt /■2tt 

d(^e""^ cos(m(/p + ipm) = yre""^'" 5mn , / dipe"'*' sin(miy9 + <^,„) = iTre""^" Sm„ , (A3) 

^0 

we find from ^ and ^ that 
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Q 



(n) 



\ Axx ^ w{x) /„ (a;) H ^ - n — ^ 



(A4) 



From ([L3), we find with the same procedure 



(n) 



2^ 



dxa;"^ w(x) 



/d ra;"+2w(a;)l ^ 



74 (x) 



(A5) 



(lAj), or that Q^"' = Q^"^ Similarly, from (111), one has 



Integrating the first term of the second integral in (|A5|) by parts and collecting terms, one immediately sees that ( A5) equals 



Q 



(n) 



r Ax {2W{x)^x-^^w{x)\ (^i^^n'^fi~f::(x)^+2n'J\xWix) (^1^ - Mfi 



(A6) 



Integrating by parts, using (|16|), and ordering terms yields that (A6) also agrees with (A4), i.e., qI"' — Q 



APPENDIX B: INFLUENCE OF MATTER INSIDE AND OUTSIDE THE ANNULUS 

We show here how the shear in the annulus i/R < \x\ < R is affected by matter inside and outside the annulus. In terms of 
the surface mass density, the shear is given by 



7(x) = - y d'yV{x^y) n{y) with V{z) = , (Bl) 

where complex notation was used, Z = zi + iz2. For a point x in the annulus, < < J?, we then define the inner and the 
outer shear by 



7in(5^) = r dyy f "^y^ , 7ou.(f) = H dyy f ^^^^ . (B2) 

^ Jo Jo - ye ^ Jr Jo - ye 

The calculation is performed by first decomposing k in Fourier modes, 



= ^ f^niy) cos(ni9 + i9„) , (B3) 



and then inserting this expression into (B2). The resulting integrals are evaluated by transforming the integration variable 
to u = e"*, di9 — — idu/it, and using the residue theorem. For the inner and outer shear, we thus obtain 

2 °° (n + l)e'''" f"^ 

7in(s) = - 7:^7772/ dy yKo(y) " / dj/ K„(y) , 



(A- 



LAJ /•OO 

7out(f)= - ^e-'''"(n-l)(A*)("-^' / dy y'^-) . 

n=2 "'^ 

(B4) 



Defining then the inner and outer multipole moments analogously by 

Qi:^ rdyy^-+'^ r A^e"'^ ^[y) = ^e"^"" /'""dy Ac„(y) , 

J Jo Jo 

Qiui := / dyy*^""' / di9 e""* k(j7) = Tre""*" / dyy^^""' , 



(B5) 



we can write the inner and the outer shear as 



ri = l 



We thus see that the shear in the annulus caused by the matter outside the annulus is fully described by the multipole 
moments of this matter distribution. The shear is afi'ected by all multipoles of the matter inside the annulus, but only affected 
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by the quadrupole and higher-order terms outside the annulus. In particular, the outer quadrupole moment causes a constant 
shear over the annulus. 



APPENDIX C: APERTURE MASSES FOR ARBITRARY APERTURE SHAPES 



We show here that the aperture mass measure ( |19[ ) can be generalized to apertures of arbitrary shape (in fact, also the 
aperture multipoles can be generalized, but they will probably be of less use). As before, we define the aperture mass as 



where w is a weight function. 



(CI) 



CI Self-similar curves 

Let the weight function be constant on self-similar concentric curves. Then define a new coordinate system by choosing a 
closed curve c(A), AG/, and setting 



y = XQ + bc{X) 



(C2) 



then the mapping between y and (6, A) is one-to-one (except at b = 0) iff cxc 7^ for all A £ /, where cxc = C1C2— C1C2. Without 
loss of generality, choose c x c > 0. The Jacobian of the coordinate transformation is J = det[d{yi,y2) /d{b, \)] = bcx c, and 
so, if depends on b only, we find from (|Cl|) that 



/•oo /• 

m{xo) = / db bw{b) (p dX (^c X K,{xo + bc{X)) 



(C3) 



We now apply the same strategy that was used for deriving Qx in the main text. Integrating (C3) by parts with respect to b 
yields 



m(f 0) = - / dbbW{b) (f) dX{cx c) ^ [xo + bc{X)) 



db 

where we have defined 

W(b) -rf db' b'w(b') . 
b Jo 

We have further assumed that 
dbbw(b) = 

to ensure that boundary terms vanish. The partial derivative dn/db can now be transformed to 

c*K _ ^ 9k, ^ ^ dhi _ ^ f dji ^ 972 j g [ 971 
db dyi dy2 \ dyi dy2 J \ dyi dy2 

1 



(C4) 



(C5) 



(C6) 



(ciC2 -I- C1C2) -—- + (C2C2 - ClCi) — — 

ob ob 



972 



+ ■ 



bcx c 



r, d-yi f 2 2\ 972 

_2C,C2— +(C,-C2) — 



(C7) 



After inserting (C7) into (C!4), the terms containing partial derivatives with respect to b and A arc integrated by parts with 
respect to the respective coordinate, making use of (C6). After collecting terms, this yields 



m{xo) = y db [2W(b) - bw(b)] j dX [- (ci£2 + C1C2) 71 + (ciCi - C2t2) 72] , (C8) 

hence the aperture mass m can be expressed directly in terms of 7. This result can be written in a more compact form by 
defining the function 



q{b) ■- 2bW{b) - bw{b) 



db' b' w{b) ~bw{b) , 



and introducing complex notation, C(A) — ci(A) -|-ic2(A). This yields 
q(bm 9 (7(j7)C*C*) 



m{xo) 



j2 

d y 



(C9) 
(CIO) 
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where C has to be taken at \{y). This generahzes the aperture mass measure as originally derived by Kaiser (1995), Kaiser 
et al. (1994) and Schneider (1996b); the generalization to non-circular apertures was first done in Squires & Kaiser (1996). 
As for circular apertures, choosing iy(6) = const, for < 6 < uB implies q{h) = for this interval, so that the mass inside the 
curve b = uB can be measured from the shear around that region. 



C2 Arbitrary curves 



We now assume that w(x) is constant on a set of nested curves of arbitrary shapes; in particular, these curves do not have to 
be self-similar. In analogy to (C2), we define a new coordinate system by 



y = xo + c{b, a) 



(CU) 



where A is a cyclic coordinate, A G J = [0, Amax], and we assume that w depends only on b. In order for (Cll) to define a new 
coordinate system, we require that the Jacobian 



j(6,A)^^\^^ 
db d\ 



(C12) 



is nonzero everywhere (except at the origin). Without loss of generality, we take J > 0. The curves will now be relabeled and 
reparameterized as follows. Let 



A{b) = - IdX&x ^ 
2 / dX 



be the area enclosed by the curve labeled with 6; then we define a new 'radial' coordinate b{b) through 
b-- 



2Aib) 



so that 

db 1 



db ^2An,ax^(6) 



dA J (b, a) 



(C13) 



(C14) 



(C15) 



We define a new cyclic parameter A (6, Aj G / by requiring that the Jacobian of the new transformation be independent of A; 
in particular, the definition of b allows to set the new Jacobian equal to b. Thus, let c{b, A) — c(b{b), A(b, A)) , then 



J(6,A) 



db dc dX dc 
db db db dX 



dX dc dfc 9A J /f f s, 
X — — = — — J(6, Aj = 6 . 



dX dX db dX 



(C16) 



where here and in the following, a prime (dot) denotes partial derivatives with respect to b (A). Equation (C16) determines 
A: using {JcU) and (ETi), we find 



A(6,A) = 



A„ 



/dA J(&(6),A) Jo 



dA J (6(6), A) 



(C17) 



These definitions of the new coordinates uniquely define at every point y the two vector fields c ' and c. We now proceed as 
in Sect. CI before, (dl) is transformed into an integral over the new coordinates (fe. A), and then integrated by parts with 
respect to b, 



m{xo) 



dbbWib) j dX^{xo + c{b,X)) 



(C18) 



with W{b) defined as in (OS), and the condition ( p6|) is applied to 10(6). The derivative of k is expressed in terms of derivatives 
of the shear, using the coordinate transform and (llCl) to obtain 



a 

db 



dX 



(C19) 



where complex notation was used again, C' = c'i+ic2. Inserting (C19) into (CIS), one obtains after another partial integration 
-(^^0) = / d^.{ - ^) 3 [(C'C)* 7] + [{C'C' - C"Cy,] } . (C20) 



We note that if one specializes c(6, A) = bk{X) to self-similar curves, (C20) reduces to (CIC) (note that (C16) implies that 
fc X fc = 1). If the weight function is chosen to be constant within a curve bok{X), one can select the family of curves such 
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that c(6A) = bk{X) for b < bo, whereas the curves are arbitrary for b > bo (except of course that the family of curves is twice 
differentiable). Then, the aperture mass (CI) can be calculated from the shear outside the curve labeled by bo- One could, 
for example, generalize the ("-statistics introduced by Kaiser (1995) to an annular region, where the inner boundary is a circle 
(e.g., situated on the cluster center), and the outer boundary is the boundary of the CCD. One can easily write down a family 
of curves which are circles for b <bo, and which smoothly deform into a square for b — B. The only difference to the standard 
(-statistics is that the vector fields c ' and c have to be calculated numerically. 
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